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Abstract 



> 

. Holonomic quantum computation is analyzed from geometrical viewpoint. We de- 

I velop an optimization scheme in which an arbitrary unitary gate is implemented 

with a small circle in a complex projective space. Exact solutions for the Hadamard, 
' CNOT and 2-qubit discrete Fourier transformation gates are explicitly constructed. 
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Recently quantum computer attracts great interests from many disciplines. It 
^ I is strongly desired to find a scheme to implement unitary gates in a physical 

I system. For this purpose it is natural to consider utilizing a quantum system 

described by a Hamiltonian H{X) that depends on external parameters {A}. 
In holonomic quantum computation (HQC) proposed by Zanardi and Rasetti 
[1,2], in contrast, the holonomy [3] associated with adiabatic change of the 
parameters along a loop in a control parameter manifold is employed to im- 
plement a unitary gate. Experimental schemes to manipulate the non-Abelian 
holonomy have been proposed [4] and their uses to realize a unitary gate are 
also proposed [5,6,7]. For efficient achievement of holonomic computation, it is 
necessary to find a loop as short as possible in the control manifold. A numeri- 
cal scheme to search the shortest loop is being developed in [8,9] to implement 
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an arbitrary gate. In this Letter we consider an ideal quantum system that 
has full isospectral parameters as control parameters. In the following we make 
exact analysis of the shortest loops that generate well-known unitary gates as 
holonomies. We would hke to emphasize that to find an exact solution has 
remained unsolved [10] even for such an idealized system. 

Here we take the relevant terminology from [11,12,13] to define our model. 
Our formulation can be compared with [14]. Suppose the Hilbert space of 
a quantum system is an A'"-dimensional complex space C"^. By isospectral 
parameters, we mean a family of unitary transformations g(X) e U{N) since 
the transformed Hamiltonians H{X) = g{X)HQg'^{\) has spectra independent 
of A. We assume that the ground state of the reference Hamiltonian Hq is 
fc-fold degenerate, taking a diagonal form 

i^o = diag(£i,£2,---,£jv) (1) 

with El = 62 = ■ ■ ■ = Ek < Ek+j {j = 1, . . . , N — k). The ground state energy 
may be put to zero, without loss of generality, to get rid of the dynamical 
phase factor. We concentrate our attention on the ground states of the Hamil- 
tonians {H{X)}. Associated with the lowest energy for each H{X), there are 
k orthonormal state vectors {]i'i(A)), . . . ,\vk{X))}. The set of k orthonormal 
state vectors is called a /c-frame and the set of all the /c-frames, 

SnA'^) = e M(7V, k; C) I V^V = 7J, (2) 

is called the Stiefel manifold. Here M(n, m; C) is the set of n x m complex 
matrices and Ik is the kxk unit matrix. The fc-frame V = {\vi{X)) , . . . , \vk{X))) 
spans a fc-dimensional subspace in C^^. The set of /c-dimensional subspaces is 
the Grassmann manifold 

GnAQ = e M(N, N; C) I P2 = P, = P, trP = k}. (3) 

The Grassmann manifold is regarded as the control manifold in the context 
of HQC. A projection map tt : S^A'^) ~^ Gn^'^) defined as 

n-.V^ 7r{V) = VVl (4) 

The group U{k) acts on SN,k{'C) from the right via matrix product as 

Sr,,k{Q X U{k) ^ SN,k{C), {V, h) ^ Vh. (5) 

Note that this action satisfies Tiiyh) = niV). Thus the set {SN,k{'C), GN,ki'C), 
Ti,U{k)) forms a principal fiber bundle with the structure group U{k). The 
group U (N) also acts on the manifolds from the left as 
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uiN) X SmAQ ^ SmAQ, (9, V) ^ gv, (6) 

U{N) X G'Ar,fe(C) ^ Giv,ik(C), {g, P) ^ ^P^t. (7) 

This action is an automorpliism of tiie principal bundle satisfying T^{gV) — 
g'K{y^g'^ . In the Stiefel manifold the connection one- form A — V'^ dV with the 
matrix elements 
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v,{X))dX'' {t,j^l,...,k), (8) 



is defined. Then the Wilczek-Zee holonomy [3] associated with the curve V{t) 
is given by U — Ve" ^ where V denotes the path-ordered product. 

We now state the main problem: Given a unitary matrix U E U{k) find the 
shortest closed loop P{t) = TT{V{t)) in the Grassmann manifold that yields the 
matrix U as its associated holonomy. This problem is called the isoholonomic 
problem and various representations of the problem have been given by Mont- 
gomery [10]. Let us consider a familiar example, a two-dimensional sphere, to 
illustrate our idea to find the solution. Suppose we parallel transport a tangent 
vector on the unit sphere along a loop starting from and ending at the North 
Pole, where the parallel transport is defined with respect to the Levi-Civita 
connection [11]. In general, the vector is get rotated from the initial direction 
when it comes back to the initial point. What is the shortest loop which im- 
plements the given rotation angle ujI The holonomy angle in this case is equal 
to the area surrounded by the loop and the problem is reduced to the so-called 
isoperimetric problem. The solution is well-known to be a small circle which 
surrounds an area of cu. We use this analogy to find a small circle solution to 
the isoholonomic problem in HQC. 

Here we describe a generalized small circle in the Grassmann manifold. Take 
/ 



Ik 




e SnAQ (9) 



as a reference point and put Pq = vr(Vo) = VqVq G ^^^^^(C) as the initial point 
of a curve in the control manifold. Taking an antihermitian matrix X G u{N), 
we define a curve V{t) = e^^Vo (0 < t < 1) in SN,k{'C) and define a curve 
P{t) = Tr{V{t)) = e*^Poe-*^ in GnA'^) by projecting V{t) into GnA'^). We 
call the curve P{t) a small circle if it satisfies -P(l) = -P(O). The connection 
evaluated along the curve V{t) is = V{tydV{t) = Vf^XVodt, from 

which we obtain the holonomy associated with the loop P{t) as 

C/ = e-/^*^ = e-^o^^°eC/(fc). (10) 
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Note that the integrand in the exponent has no ^-dependence and Jq dt simply 
yields a factor unity. 

Suppose that we are to implement a unitary gate Ug^te by HQC. Our task is 
to find a control matrix X e u(A^) such that the holonomy (10) reproduces 

f/gatc- In general, computing the holonomy U for a given X is easy but the 
inverse problem is considerably difficult; we need to find X for a given U while 
keeping the closed loop condition P(l) = -P(O) satisfied. Here we report some 
of the exact solutions of this inverse problem for several important gates. 

Define the logarithmic matrix fl £ u(/c) of the gate such that f/gatc = If the 
eigenvalues of f/gatc are (e"*'^^, e~*'^^, . . . , e""^*), fl has eigenvalues {iui,iu!2, ■ ■ ■ , 
icuk) in the range —n<ujj<Tr and the same eigenvectors as f/gate- Then the 
problem (10) reduces to find X e u(A'") such that 

v^xvo - n. (11) 

The general solution to this problem takes the form 

X = I (12) 



-W^ z 



with matrices Z e u{N - k) and W e M{k, N - k;C). Montgomery [10] has 
shown that any optimal solution necessarily satisfies Z = 0. The remaining 
problem is to find the matrix W that satisfies the loop condition P(l) = -P(O). 
This condition demands that be of the form 

e-=\*°\. (13) 
** 



where * and ** are nonvanishing matrices. The solution W = is not ac- 
ceptable since it gives a constant curve P{t) = Pq, leaving the Hamiltonian 
unchanged. To seek a nontrivial solution we introduce a penalty function p{X), 
which measures the norm of the off-diagonal-block elements of e^, by 



N 



=1 j=k+l 



(14) 



being the set of orthonormal basis vectors of C^. Then the closed loop 
condition (13) is rephrased as an equivalent condition ]9(X) = 0. It is apparent 
that p{X) > by definition and this nonnegativity is suitable for numerical 
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search of the zeroes of p{X). Through numerical studies we have arrived at a 
method to construct systematically exact solutions of the equation p{X) — 0. 

Here we describe the method briefly. We restrict ourselves to the cases such 
that = A; + 1 for simplicity. In this case the Grassmann manifold GN,ki'C) 
reduces to the complex projective space CP^~^ = Gn,n~i{'C)- A given unitary 
gate Legate = has a set of eigenvectors and eigenvalues {{uj, e"""^^ ) | uj G C'^, 
ujj e M, —n < uij < n, fluj = iUjUj, {uj,ui) = 5ji}. Then choose a pair of 
eigenvector and eigenvalue, {u^,u!fj), and substitute 

W(^n,n) = a(^n,n)Uf„ a(^f,,n) = (Sttti + uj^) {27rn - Uf,) (15) 



into of X in (12) for n = ±1, ±2, . . .. Then we obtain the solution 
that satisfies p(X(^^„)) — 0. Therefore, the integer n generates a family of 
solutions and there are inequivalent families of solutions as many as different 
eigenvalues of f/gatc- The proof of the above method is lengthy and will be 
published elsewhere. In the rest of this Letter we show results exemplifying 
this method. 



We first consider the Hadamard gate 

(16) 




which has two eigenvalues and corresponding eigenvectors 



1, Ui 



cos g 
sin ? 



-1, U2 



— sm 
cos ? 



J 



(17) 



Hence there are two families of exact solutions of the equation p{X) — 0. 
We construct the first family by substituting the eigenvector and eigenvalue 
{ui,uji) into (15). Then the control matrix X of the first family is 



/ 



(1) 

H 



ITT 



sin^? 



— sin ? cos ? 



ne 



cos ■ 



- sm I cos I ne'" cos | 
cos^ I ne*^ sin | 
ne~*^ sin 5 



(18) 



where 6 is an arbitrary real number, which parametrizes unitarily equiva- 
lent solutions. The integer n counts the winding number of the loop P{t) ~ 
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e*^Poe *^ as seen below. The eigenvalues and eigenvectors of X^^ are easily 
found to be 



Xi — ITT, \Xi) — 



( ■ A 

- sm f 



cos 




e cos ^ 



X2 = ^n7^, \X2) 



V2 



e'^ sin f 



(19) 



xs^-mn, 1x3) 



V2 



( e'^ cos ? ^ 



e sm ^ 
-1 



Using the spectral decomposition of we can calculate its exponentiation 
e*^, from which we obtain the penalty function as 



p{tx'^^)^sin^{n'Kt). 



(20) 



This clearly shows that the loop P{t) passes through the initial point Pq 
\n\ — 1 times in the interval < t < 1. Accordingly, the number n describes 
how many times the loop winds as t changes from to 1. It is clear that 
n = ±1 is nontrivial and optimal in this family of solutions. 

The second family of solutions is constructed by substituting the second eigen- 
vector and eigenvalue (1*2, 1^2) of C/r into (15) as 



(2) 



sin^? 



-smfcosf 
-a„e~*^ sin f a„e"*^ cos f 



sin f cos f — OnC*^ sin f ^ 



cos^ ? 



(21) 



where a„ = \/4n2 — 1/2. The eigenvalues and corresponding eigenvectors of 
X'i^ are 



6 



xi = 0, \xi) = 



( A 

COS I 
sin ^ 



x^ — tTT I —n + 



\ " / 



X2^i7r(n+^^, \x2) 



—bne sin ^ 

6„e^^ cos? 



(22) 



Cy^G- COS g 



where 6„ = ^^(277- + l)/4n and = y'' (2n — l)/4n. This time the penalty 



function is evaluated as 



(23) 



The integer n is again interpreted as the winding number of the loop and the 
choice n = ±1 provides the optimal solution in this family. 

We have shown that there are two families of exact solutions to the holonomic 
implementation of the Hadamard gate. To determine the shortest loop we 
need to calculate length of each loop employing the Fubini-Study metric of 
the Grassmann manifold. Consequently, the velocity of the point P(t) and 
therefore the length of the loop are proportional to the norm of the 

matrix W . These norms are evaluated to be 

\\WS\\ = 7r|n|, \\W^\\ = ^V4n2 - 1. (24) 



Therefore, we conclude that the simple loop solution n 
family is in fact optimal in the whole class of solutions. 

Next we turn to the CNOT gate 

( \ {) ^ A 

10 

1 

10 



±1 in the second 



CNOT 



(25) 



7 
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For this case there are two famihes of solutions. The first one is given by 



CNOT — ™ 





























1 nd2e 








_ 1 

2 

1 

2 







1 

2 



72 



nd^e 



ndie'^' 
-^ndse'^^ 



J 



with real numbers {d\^d2^d^ with a constraint o?^ + (i| + (i| = 1. The param- 
eters {^1,^2,^3} are arbitrary real numbers. The integer n again counts the 
winding number of the loop. The second one is 



(2) 

CNOT 





















1 

2 



-^OnC ''^ Ts^ne' 








(27) 



with a„ = \J ^v? — 1/2. By similar comparison of the norms ||VFcnot|| &s in 
(24), we conclude that the choice n = ±1 in the second family yields the 
optimal loop. 

Our final example is the 2-qubit discrete Fourier transform (DFT2) gate 



DFT2 



( \ \ 1 1 ^ 



1 i -1 -i 
1-11-1 
1 -i -1 i 



(28) 



There are three famihes of the solutions for this gate. The first family takes 
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the form 



(1) 



DFT2 



i nwi 
nw2 
-J riW3 



—J — I nw4 



^^nw^ ^ 



(29) 



with the parameters 



= W4 



(30) 



with a constraint /i + /I = 1- The second family is 



(2) 

DFT2 



- 



1 

4 

_ 1 

2 

_ 1 
4 









oS'"^' 








(31) 



with = y/lGn"^ — 1/4. The third family is 



/ 



X 



(3) 

DFT2 



_1 

4 

1 
4 

1 
4 

1 
4 

" 2 



2 



2^n^ 2^n^ 



2 ^TI.^ 



a„e 
One 



(32) 



with Un — -\/4n^ — 1/2. The norm of the matrix W with n = ±1 for each 
family is now evaluated as 



I|w^St2| 



I|W^St2| 



I|w^St2| 



(33) 
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Thus the simple loop (n = ±1) in the third family gives the optimal control. 

In this Letter we considered an ideal system which has full isospectral param- 
eters as control parameters. We found exact implementation of Hadamard, 
CNOT and DFT2 gates with a small circle in the complex projective space 
([^pN-i _ Implementation of larger-qubit gates is under progress 

and will be published elsewhere. 

A reahstic system has a restricted control manifold M and a control map 
f : M ^ GN,k{'C). For physical realization of HQC it is required to find an 
optimal loop for the holonomy in the pullback bundle f*{SN,k{'C)) as discussed 
in [2,4,5,6,7]. The exact solutions we have obtained here are pulled back by /* 
to be loops in M. Detailed analysis of physical realization is beyond the scope 
of this Letter and will be pubhshed elsewhere. 
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